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Neutrino interactions with lead produce neutrons in numbers that depend on neutrino energy and 
type. A detector based on lead perchlorate, for example, would be able to measure the energy 
deposited by electrons and gammas in coincidence with the number of neutrons produced. Sorting 
the electron energy spectra by the number of coincident neutrons permits the identification of the 
neutrino type that induced the reaction. This separation allows an analysis which can determine the 
temperatures of z/ e and v e from a supernova in one experiment. The neutrino reaction signatures of 
lead perchlorate and the fundamentals of using this material as a neutrino detector are described. 

PACS numbers: 97.60.Bw, 25.30.Pt 

I. INTRODUCTION 



Recently a number of groups have expressed interest 
in using Pb as a target for neutrino interactions to study 
supernovae or oscillations ^j, and this inspired an 
estimate of the cross section [Q. The interest arises be- 
cause of the large cross section and the low relative cost 
of Pb. As a result, additional cross section calculations 
were done recently by Fuller, Haxton, and McLaughlin |5|] 
(hereafter referred to as FHM) and Kolbe and Langanke 
HD (hereafter referred to as KL). 

The interesting neutrino interactions with Pb consist 
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Fig. [l] shows the energetics of these transitions. The 
number of neutrons emitted (0, 1, or 2) depends on the 
neutrino energy and whether the transition is induced by 
a charged current (CC) or neutral current (NC) interac- 
tion. The nuclear physics of this system is described in 
detail in Rcfs. |.|. 

An ideal lead-based neutrino detector would have an 
appreciable density of lead atoms and the capability of 
detecting the electrons, gammas and neutrons produced 
in the reaction. Lead perchlorate (Pb(C104) 2 ) has a very 
high solubility in water (500 g Pb(C10 4 ) 2 / 100 g H 2 §) 
and the mixture is transparent. This transparency raises 
the hope that a Cerenkov detector can be assembled. 
Additionally, the presence of a neutron moderator (hy- 
drogen) and a neutron capture nucleus ( 35 C1) provides 
a method for observing the neutrons. The high-energy 
(8.6 MeV), neutron-capture 7 rays from CI would Comp- 
ton scatter in the fluid and be observed via the Cerenkov 
light of the recoil electrons. 

Supernovas emit on the order 10 58 v of all types. The 
average energies of neutrinos emitted by a supernova fol- 
low a hierarchy: E Vt < E^^ < E v T where v^. T indicates 
either /i or r neutrinos and their anti-particles. The pre- 
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FIG. 1. The level scheme of 208 Pb - 208 Bi system. The lev- 
els labeled with GT indicating the Gamow- Teller resonance, 
IAS indicates the isobaric analog state, and 1st indicating the 
states populated by first forbidden transitions. The 1, 2, and 
3 neutron emission thresholds are indicated by the labels In, 
2n, and 3n respectively 



dieted values of the average energies for the three neu- 
trino classes differ a great deal and are estimated to be 
E v = 11 MeV, E P = 16 MeV, and E v = 25 MeV 
||. 

Although not strictly thermal, the supernova neutrino 
energy spectra can be described well by a thermal distri- 
bution, 



fu 
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where /„ is the normalized neutrino spectral shape, T„ is 
an effective neutrino temperature, E v is the neutrino en- 
ergy, and 77 is the degeneracy parameter (chemical poten- 
tial divided by T v ). For v e and v e , the supernova model 
predictions for the flux is well fit with 77 = 3, whereas 



for 



the flux is better described with -q = 0. The 



normalization factor F 2 (?y) has the values F 2 (0) = 1.803 
and F 2 (3) = 18.969. The total neutrino fluence F v from 
a supernova can be estimated by 
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TABLE I. The neutron capture cross section for the var- 
ious species in Pb(C104)2 the number of species in solution 
relative to Pb for an 80% solution concentration and the neu- 
tron capture rates relative to hydrogen. 
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where Eb is the released energy and D is the distance to 
the supernova. 

The supernova v e spectrum is too soft to produce 
a large number of multiple-neutron inverse-beta-decay 
events in Pb. However, if the higher energy v^. T oscil- 
late into v e , the hardening of the spectrum will greatly 
alter the expected number of 2-neutron events . Thus 
the response of a Pb detector is intriguing because if T, yc 
is found to be larger than T Pe , it would be strong evi- 
dence that neutrino oscillations are taking place. How- 
ever, the uncertainties in the neutrino energy distribution 
as predicted by the supernova models may complicate the 
interpretation of a measurement of the spectral parame- 
ters. 

The FHM paper demonstrated that the NC and CC 
interaction rates in Pb are sensitive probes of T„ „ r , 
T„ c , and T Pc . This paper further discusses the details of 
studying supernova neutrinos with Pb. In particular the 
use of the added data arising from observing the product 
electrons and gammas in coincidence with neutrons is de- 
scribed. These data permit the separation of the NC and 
CC rates and a further division of CC v e and v e events. 
This collection of data can be interpreted in terms of the 
3 temperatures of interest. The CC event division is dis- 
cussed here in detail. Although any Pb-based experiment 
that measures electrons, gammas and neutrons might ap- 
ply the techniques described herein, the example of using 
lead perchlorate as a target material is discussed. 

II. PHYSICAL AND OPTICAL PROPERTIES 

To build a reasonably large detector viewed by photo- 
multiplier tubes from the periphery, the attenuation of 
the Cerenkov light must be minimal. Therefore the 



transmission of light through the medium is critical. To 
understand the transmission, measurements were made 
at several wavelengths through several concentrations of 
Pb(C104)2 pxj ]- The maximum achievable transmission 
is still under study with attentuation lengths exceeding 
2.5 m having been measured in an 80% solution. A large 
mass detector may require segmentation to minimize ad- 
verse effects from light loss. 

The Cerenkov response of the medium depends on the 
index of refraction, which for Pb(C104)2 is 1.5 for an 
80% solution by weight. The density of this solution is 
2.7 g/cm 3 , the corresponding 208 Pb number density is 
1.7 x 10 21 /cm 3 and the hydrogen number density is 3.6 x 
10 22 /cm 3 . Some additional pertinent data for Pb(C104)2 
are given in Table 

The number of Cerenkov photons emitted by an en- 
ergetic electron in Pb(C104)2 is about 185/cm |ll[] in 
the wavelength region of interest for phototubes. The 
stopping power of 80% (50%) Pb(C10 4 ) 2 is about 0.2 
cm/MeV (0.33 cm/MeV). Hence a 15 MeV electron will 
emit a total of about 550 (920) photons in 80% (50%) 
Pb(C104)2- Due to the modest light yield, the resolution 
(«20%) will be modest. The resonance widths will also 
contribute an uncertainty in relating the electron energy 
to the incoming neutrino energy by a comparable amount 
(few MeV). 

Because solutions of Pb(C104)2 contain a great deal 
of water and CI, the neutron thermalization and capture 
time is of the order of 10-100 /is. This time scale is short 
compared to the v detection rate during a supernova, 
but long enough to identify whether there are 0, 1, or 2 
neutrons in coincidence with a primary e~ or 7. 

III. CROSS SECTION MODELS 

The FHM Q paper pointed out the importance of the 
208 Pb - v cross section and the production of neutrons. 
They noted that the total number of neutrons produced 
by neutrino interactions with Pb is sensitive to the ef- 
fective temperature of the supernova v e energy spectrum 
especially for reactions which produce multiple neutrons. 
However, they did not consider detection schemes that 
could measure the energy released in the form of elec- 
trons and gamma rays (hereafter referred to as the elec- 
tromagnetic energy). Neutral current events release a 
large number of neutrons but very little electromagnetic 
energy. Hence NC events can be isolated from CC events 
by measuring the electromagnetic energy in coincidence 
with the neutrons. 

The papers by FHM and KL provided effective cross 
sections for neutrino-lead interactions averaged over the 
neutrino energy spectra and summed over all the prod- 
uct nucleus states. It is just these distributions, however, 
which are needed to simulate the response of a detec- 
tor and to optimize its design parameters. This section 
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TABLE II. A description of the levels used in the simula- 
tion. The 2 neutron rate from the NC 1st forbidden transition 
is not critical to the analysis. The neutron number branching 
ratios used here were motivated by results from FHM. 
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0.45 





NC first forb. 


14 MeV 





~ 1 


~ 



I ■ I 



1000 - 



> 





o 



1100 - 



10 - 



LU 

is 

T5 



1 - 



0.1 - 



I ■ I 




T 



T 



T 



- 80x10 



- 60 «-* 



> 

cr 

40 J 



- 20 



J— 



40 80 120 
Energy (MeV) 

FIG. 2. The CC cross section as a function of energy plot- 
ted on the left axis (solid line). A normalized thermal neu- 
trino energy spectrum (T„ c = 2.76 MeV, r\ = 3.) is plotted 
on the right axis (dashed line). Note the cross section scale is 
logarithmic whereas the spectrum scale is linear. 



summarizes the nuclear physics model implemented to 
calculate these distributions. In the following section, 
these distributions and their use in analyzing the neu- 
trino spectra are described. 

Specifically, the energy dependence of the five primary 
transitions of interest as identified by FHM is calculated: 
two for the NC interaction and three for the CC inter- 
action. The parameters describing these transitions are 
given in Table 0. For the isobaric analog state (IAS) 
CC and Gamow-Teller (GT) NC and CC transitions, the 
analytic expressions and matrix element values from Ref. 
H are used. For the first forbidden transitions, the cross 
section can not be written in closed form. Software to 
calculate that contribution to the cross section as used 
by FHM was obtained Q] and hence the cross section 
calculations presented here agree with FHM by design. 
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FIG. 3. The CC cross section as a function of energy plot- 
ted on the left axis (solid line). A normalized thermal neu- 
trino energy spectrum (T„ c = 6.27 MeV, r] — 3.) is plotted 
on the right axis (dashed line). Note the cross section scale is 
logarithmic whereas the spectrum scale is linear. 



Figure 2 shows the cross section overlaid on a nor- 
malized thermal flux distribution (X'„ c =2.76 MeV, r?=3). 
Figure 3 shows the same for a higher neutrino tempera- 
ture (2^=6.27 MeV, T}=3). The origin of the large in- 
crease in the interaction rate with temperature is clear 
from a comparison of these plots. The onset of the 
Pb cross section lies between these two temperature ex- 
tremes. As a result, the interaction rate in Pb is highly 
sensitive to T„ e . This is especially true for the first forbid- 
den transition which leads to the emission of 2 neutrons. 



Tables III and IV show some cross sections and rel- 
ative interaction rates for the various components of 
Pb(C104)2- It is clear that the interactions with Pb 
and H dominate. Contributions from CI, O, and elas- 
tic scattering of electrons are at the few percent level. 
The CC interaction of v e on protons is a major contribu- 
tion because the hydrogen density is large in Pb(C104)2 
solutions. For other isotopes found in Pb(C104)2 but 
not listed in the table, FHM note that the cross section 
models vary like N-Z or NZ/A. For Pb, these dependen- 
cies are small for the naturally occurring isotopes. Hence 
one can anticipate that their response to neutrinos will 
be similar. This is not the case for CI. However, the 
difference in cross section between 35 CI and 37 C1 is still 
probably smaller than other uncertainties and thus the 
CI interaction rate is small. 
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TABLE III. Cross sections for v e interactions on several 
isotopes, for various energy spectra and authors. The column 
labeled by DAR indicates an average of the cross section over 
the muon decay-at-rest energy spectrum for electron neutri- 
nos. The column labeled SN is for a Fermi-Dirac distribution 
of energies with 77 = 3 and T v = 6.266 MeV (E v ~ 25 MeV). 



The cross sections are in units of 10 
various other cross sections are given. 



cm . For reference 



Isotope 



' 2us Pb(!/ e ,e-)X 

7 Cl(!;e:0 37 Ar 

35 ci(i> c ,ox 

35 Cl(P e ,e + )X 
16 0(^ c .e-)X 
16 OK,<)X 
D(j/ e ,e~)X 
H(p e ,e+)n 
^ e (e,e')i/ e 




DA 



59.6 



0.53 








a This value was calculated by the author using the formalism 
of Ref. §. 

b This value is for 77 = and T„ = 8.0 MeV. 

TABLE IV. The neutrino interaction rates for 30 MeV v e 
for the dominant reactions. The values are given relative to 
208 Pb. The cross section value for Pb corresponds to that of 
FHM. 



Element 


Relative Rate 


2m ~Pb{v e ,e) 
37 Cl(^,e) 
ls O(i-„c) 
H(p e ,e+)n 
f e (e,e')i/ 


100 
2.0 
0.6 
32 
3.1 



IV. THE ELECTROMAGNETIC ENERGY 

A detector based on lead perchlorate would be able to 
measure the energy deposited by e~ and 7 in coincidence 
with the number of neutrons produced. The energy de- 
posited in the detector due to e" and 7 is here referred 
to as the electromagnetic energy to contrast it with that 
resulting from neutron capture. To identify a particular 
class of neutrinos, one sorts the electromagnetic energy 
spectra by the number of coincident neutrons. CC events 
always produce an electron which can have substantial 
energy in coincidence with neutrons. NC events produce 
either a high-energy 7 with no neutrons or neutrons with 
little or no 7 energy. 

First consider the separation of NC from CC events. 
The product nucleus from the neutrino interaction can 
produce 0, 1, or 2 neutrons. (Particle decays other than 
neutrons (alphas and protons) of the product nucleus 
states have a much lower branching ratio (see KL) and are 
ignored here.) When the interaction does produce a neu- 
tron, the neutrons carry away the available energy and 
there is little electromagnetic energy produced. Hence 
for the NC events, the electromagnetic energy is only 
appreciable when neutrons are emitted. Thus events 
which do have appreciable electromagnetic energy in co- 
incidence with neutrons arc indicative of 208 Pb - v e or H 



- P e CC interactions. 

Next consider the separation of the CC events into v e 
or v e events. The 2-neutron event rate is dominated by 
ivinduced CC transitions to the first forbidden level in 
Bi. Since the electromagnetic energy for the CC reactions 
is related to the incoming neutrino energy, the 2-neutron 
event spectrum provides data on the incident v e energy 
spectrum. In contrast, the IAS and GT transitions in- 
duced by v e produce mostly 1-neutron events. Further- 
more, although D e do not have an appreciable CC cross 
section on Pb, they have a large cross section on H. The 

1- neutron spectrum is therefore generated by reactions 
induced by v e and v e . But those are the only reactions 
that contribute significantly to the 1-neutron electromag- 
netic energy spectrum and therefore the 1-neutron spec- 
trum provides data on the i> e energy distribution. The 

2- neutron and 1-neutron spectra thus relate directly to 
the spectral features of the incident neutrino flux. Figs. 
|] through ^ show spectra of the electromagnetic energy 
for events in coincidence with 1 or 2 neutrons for two dif- 
ferent neutrino energy distributions. These spectra have 
been convolved with an estimate for the detector resolu- 
tion. 

For NC events, only GT excitations can produce an 
event with no neutrons as the resonance straddles the 
neutron separation energy. A 7 ray of approximately 7.3 
MeV is then emitted. This would produce a line in the 
0-neutron electromagnetic energy spectrum. NC tran- 
sitions can also produce 1-neutron events and a mod- 
est number of 2-neutron events. As there will be lit- 
tle electromagnetic energy in coincidence, the 1-neutron 
events will be interpreted as an 8.6-MeV line in the 0- 
neutron spectrum. The strengths of these two lines in 
the 0-neutron spectrum provide a measure of the NC 
rate. Due to the resolution, these two lines will most 
likely be blended. These rates will be large even for low- 
temperature profiles. 

One background to the CC 1-neutron spectrum comes 
from 2-neutron NC events. One of the 2 neutrons would 
be taken as a neutron but the second would be mistaken 
for electromagnetic energy. For T„ ,„ T =6.27, T Vc =2.76, 
?7=3, the NC rate of 2-neutrons events will be comparable 
to the CC rate because, although the NC cross sec- 
tion is much smaller, there are 4 times as many v^ T as v e 
and their temperature is much higher. Thus at low Ti, E , 
this background must be separated. Note that because 
there are several 7 sharing the energy and photons Comp- 
ton scatter more than once before losing all their energy, 
the Cerenkov light produced by neutrons capturing on CI 
will be isotropic in comparison to electrons from inverse 
beta decay. This could lead to a technique for separating 
this background. Otherwise an electromagnetic energy 
threshold of 10 MeV would eliminate it. 

With the spectra in hand, one deduces the v temper- 
ature by relating it to the average electromagnetic en- 
ergy from the 1-neutron spectrum (E\ ) and the 2-neutron 
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FIG. 4. The energy observed in the detector due to elec- 
trons from the CC interaction that are found in coincidence 
with 2 neutrons (E 2 ). This plot is for T„ e = 2.76 MeV, r\ = 
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FIG. 6. The energy observed in the detector due to elec- 
trons from the CC interaction that are found in coincidence 
with 1 neutron (Ei). This plot is for T„ e = 2.76 MeV, T De = 
4.01 MeV, r\ — 3. The dashed line represents the contribution 
due to i> e on H. 
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FIG. 5. The energy observed in the detector due to elec- 
trons from the CC interaction that are found in coincidence 
with 2 neutrons (E 2 ). This plot is for T„ e = 6.27 MeV, 77 

:s. 
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FIG. 7. The energy observed in the detector due to elec- 
trons from the CC interaction that are found in coincidence 
with 1 neutron (Ei). This plot is for T„ c = 6.27 MeV, Tp e = 
4.01 MeV, 77 = 3. The dashed line represents the contribution 
due to i> e on H. 
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Electron Anti-Neutrino Temperature (MeV) 

FIG. 8. A contour plot of E\ (solid lines) and E2 (dashed 
lines) as a function of the neutrino temperatures. In both 
cases, the average energy was calculated for events depositing 
at least 5 MeV. 

spectrum (E 2 ) (Fig. ||). One also counts the total num- 
ber of 1-neutron events (Ni) and 2-neutron events (-/V2) 
in these spectra (Fig. ||). The ratio (R=N 2 /Ni) depends 
strongly on the temperatures of the two types of neutri- 
nos but is insensitive to uncertainties in the supernova 
distance scale. Using these parameters (E~i, E 2 , N 2 , N\, 
R), the pair of neutrino temperatures that fit the data 
are determined. 

If the electron neutrino temperature is clearly greater 
than that of the electron anti-neutrino temperature, it 
would be considered strong evidence for v e — or v e — v T 
oscillations. This condition is satisfied by the top left 
hand section of these two-temperature plots. It should 
be noted that the contour plots were done for thermal 
spectra. If v oscillate, the spectral shape could be com- 
plicated depending on the oscillation conditions and this 
analysis would have to be extended. 

To use this set of figures to determine the temperatures 
with some precision, one must estimate the uncertainty in 
the measured parameters. These uncertainties are mostly 
determined by the number of events observed. Regions 
in Figs. H through [To] are defined by the values of these 
parameters and their uncertainties. The boundaries of 
these regions indicated on the axes provide the corre- 
sponding temperature ranges that describe the data. For 
a commonly considered supernova at 10 kpc with 3 x 10 53 
ergs of released energy and a 10 kt-Pb detector (24 kt of 
lead perchlorate) , one expects to measure E 2 = H-0 ± 
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Electron Anti-Neutrino Temperature (MeV) 

FIG. 9. This is a contour plot of the count rates as a func- 
tion of the neutrino temperature for a 24 kt Pb(C104)2 de- 
tector. The solid lines are for the number of 1-neutron CC 
events and the dashed lines are for the number of 2-neutron 
CC events. A threshold of 5 MeV was used. 




Electron Anti-Neutrino Temperature (MeV) 

FIG. 10. This a contour plot of R as a function of the 
neutrino temperatures. R is the ratio of two-neutron events 
to one-neutron events. A threshold of 5 MeV was used. 
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0.5 MeV, E 1 = 19.6 ± 0.2 MeV, N x = 1530, N 2 = 82, 
and R = 0.053 ± 0.006 for non oscillating spectra {T Vc 
= 2.76, 7] = 3, E Ve = 11; T Pe = 4.01, r) = 3, S Pc = 16)". 
The uncertainties for the average energies were estimated 
by calculating the variance of the distribution and divid- 
ing by the square root of the number of observed counts. 
The resulting statistical uncertainty in the temperatures 
is small and the precision improves if T„ e is larger because 
the event rate is higher. 

V. SUMMARY 



TABLE V. Some cross section comparisons between FHM 
and KL. The cross sections are given in units of 10 -40 cm 2 . 
In most cases the two calculations have fair agreement. The 
selection here was chosen to emphasize the differences. 



Transition and 


FHM 


KL 


spectral input 






Pb{u e ,e)X T„ e =4.0 MeV, n=3. 


6.9 


6.7 


Pb(i/ e ,e)X T„ e =8.0 MeV, 7j=0. 


58 


43 


Pb(i/,i/')X T„ e =4.0 MeV, r?=3. 


0.66 


0.23 


Pb(v,v')X T„ e =8.0 MeV, rj=Q. 


4.5 


1.4 


Pb(i/,e)X tt^ decay at rest 


91 a 


59.6 



a This value was calculated by the author using the formalism 
of FHM. 



The Sudbury Neutrino Observatory (SNO) |19| and 
Super-Kamiokande |2fj|] experiments are operating. In 
the event of a supernova, Super-Kamiokande will ob- 
serve a large number of v e -K events (~ 10000 for Tp e 
= 4.01, rj = 3) and SNO will see a number of v e -T) events 
(~ 90 for T Ue — 2.76, r\ = 3). The average energy mea- 
sured for these two reactions could be compared in a 
two-temperature plot in a similar analysis as described 
above. Due to the smaller number of events, the uncer- 
tainty comparing these energy distributions would be de- 
termined by the SNO measurement. For these v tempera- 
tures, a 24 kt Pb(C104)2 experiment would count a num- 
ber of 2-neutron v e events comparable to the number of 
SNO D(^ e ,e~)pp events. Thus the average energy anal- 
ysis of a Pb(C104)2 experiment would have a sensitivity 
similar to the SNO-SK comparison. But the total num- 
ber of counts in a Pb(C104)2 experiment is much more 
sensitive to the neutrino temperatures. SNO (Pb(C104)2 
) would observe about 3 (60) times as many D(i/ e ,e _ )pp 
(2-neutron) events at T Ve = 6 than at 2.76 MeV. This 
extra handle on the data adds a powerful redundancy to 
the analysis. 

It is obvious from Table that there are still large 
uncertainties in the nuclear physics of these interactions. 
Furthermore, reactions on the other Pb isotopes have not 
been considered although they would tend to increase 
the signal. Therefore the response of such a detector 
should be measured in a prototype. The v e spectrum 
from stopped decay is comparable to that from su- 
pernovae. So the 2-neutron and 1-neutron spectra could 
be measured at a beam stop facility pl[| without the in- 
fluence of v e . Since the £/ e -proton reaction is well under- 
stood, it does not need additional study and the measure- 
ment would focus on what is currently poorly known: the 
nuclear physics parameters of the lead reactions. Thus 
the theoretical uncertainties in the cross section could be 
diminished by a beam stop experiment. 

The details of the analysis presented in this paper will 
alter as the nuclear physics of this system becomes bet- 
ter known. But the general conclusions generic to lead 
detectors will not. It is very powerful to measure not just 
the total number of neutrons produced but the 2-neutron 



and 1-neutron rates separately. Furthermore the obser- 
vation of electrons and gamma rays in coincidence with 
the neutrons can permit the separation of NC and CC 
interactions. If one observes a supernova with a Pb de- 
tector, these combinations of data can be analyzed to 
estimate the temperatures of the neutrino spectra. 
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